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ABSTRACT 

A  method  for  analyzing  in  the  frequency  domain  the  performance 
of  linear  networks  containing  nonlinear  resistors.  This  method 
is  applied  to  the  evaluation  of  the  frequency  performance  of  a  re- 
actively  terminated  mixer. 
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FREQUENCY  DOMAIN  ANALYSIS 
OF  A  CLASS  OF  NONLINEAR  NETWORKS 


The  evaluation  of  the  frequency  performance  of  networks  containing  nonlinear  devices  is 
normally  very  involved.  This  report  describes  a  procedure  by  which  the  frequency  domain  per¬ 
formance  of  networks  containing  devices  that  may  be  represented  as  nonlinear  resistors  can  be 
evaluated  with  any  desired  accuracy.  Since  only  general  statements  can  be  made  abnui  tin*  appli¬ 
cation  of  this  method  to  the  whole  class  of  networks  considered,  the  detailed  analysis  of  a  bal¬ 
anced  mixer,  together-  with  its  frequency  selective  terminations,  is  presented  as  an  example. 

The  Problem 

The  general  problem  considered  is  a  network  containing  nonlinear  resistors  controlled  by 
one  or  more  independent  voltages  and/or  currents  embedded  in  a  linear*  frequency  invariant  net¬ 
work  (Fig.  1).  Also,  ports  1  and  Z  are  assumed  to  be  the  input  and  output  ports  of  a  2 -port  X 
(dashed  lines,  big.  1).  Later,  the  problem  of  embedding  N  in  a  general  linear  network  including 
reactive  elements  is  considered. 

h-  s-imi| 


I  I 
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Fig.  1.  Frequency  invariant  network  containing  nonlinear  resistors. 

Some  Theor  etical  C'onsiderat  ions 

‘The*  nonlinear  resistors  are  assumed  to  be  monotonic  and  representable,  over*  the  frequency 
band  of  interest,  by  some  curve  or-  family  of  curves  (Fig.  2).  The  single  curve  type  [i.e.,  Lig.  2(al| 
will  bo  a  two-terminal  device  (e.g.,  a  diode).  If  the  device  possesses  three  terminals,  with  one 
used  in  the  fashion  of  a  grid  on  a  triode  or  gate  on  an  I  ET,  the  multiple  curve  type  description 
[Fig.  2(b)]  is  necessary.  By  analogy  this  scheme  can  be  extended  to  cover  more  control  param¬ 
eter. s.  The  currents  flowing  through  the  nonlinear  resistors  will  be  referred  to  as  f.(v),  f.(v v^), 
etc.,  respectively  (Fig.  2)  to  distinguish  them  from  currents  in  the  linear-  porlion  of  the  network. 

\ssume  for  the  time  being  that  only  m  two-terminal  nonlinear  resistors  a  re  preseni  in  die 
network.  If  these  resistors  are  extracted  to  form  m  ports  and  if  the  n  input  ports  are  supplied 
wilh  voltage  sources  Y  .  fed  through  resistors  Rfi.,  Fig.  3  results.  It  will  be  assumed  for  the 

time  being  that  the  \  ^  are  DC'  voltages.  The  reason  for  this  will  become  apparent  later.  NR 
st  th 

stands  for  nonlinear  resistance  i.  The  current  at  the  i  output  port  is  then  f.(v.)  where  is 
the  voltage  across  the  port.  Defining  a  2n  x  2m  transmission  matrix  between  the  n  input  and 
m  nonlinear  resistor  ports  with  elements  A,  B,  C,  and  D  (each  of  which  is  an  n  x  m  matrix  con¬ 
taining  constants  only). 
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Fig.  2.  Typical  nonlinear  resistance  curves. 


Fig.  3.  (n  +  m)-port  linear  network 
terminated  in  nonlinear  resistors  at 
the  m -ports. 
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Vs  ( R g C  +  A)  v  +  (R  D  +  B)  f  (v) 


(1) 


where  V  is  an  n  X  1  column  vector,  v  and  £ (v)  are  m  x  1  column  veetors  and  R  is  an  n  x  n 
diagonal  matrix  with  resistor  Rg.  in  the  i  X  i  position  (for  i  =  1,2...  n). 

If  three-terminal  nonlinear  devices  are  now  considered  as  well,  then  for  each  such  device 
wherever  there  was  one  output  port  in  Fig.  3  there  will  now  be  two.  This  scheme  is  illustrated 
in  Fig.  4,  It  is  assumed  that  no  current  flows  into  the  control  port.  This  can  obviously  be  done 
with  no  loss  in  generality  hy  merely  assuming  the  resistance  (if  any)  across  this  port  as  being 


Fig.  4.  Transformation  of  a  single 
to  a  double  controlled  nonlinear 
resistor. 
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absorbed  into  the  main  network.  So  if  there  are  k  of  these  th ree -terminal  devices  present,  the 
number  of  output  ports  will  be  increased  to  (m  +  k).  An  equation  essentially  in  the  same  form  as 
that  given  in  Fq.  (1)  will  still  hold;  however,  the  transmission  matrix  elements  will  be  n  x  (m  +  k) 
matrices  and  £(v)  will  he  a  function  of  two  voltage  variables,  v.  and  V(  .,  for  each  of  the  three- 
terminal  devices.  If  devices  with  more  than  two  control  parameters  are  used,  the  foregoing 
principles  still  hold  merely  increasing  the  order  of  the  transmission  matrices  and  appropriately 
('hanging  £(v).  Also,  the  controlling  parameters  need  not  be  voltages,  they  may  he  any  mixture 
of  currents  and  voltages.  No  matter  what  the  complexity  of  the  devices,  however,  an  equal  ion 
of  the  form  shown  in  Fq.  (1)  will  always  be  obtainable  by  standard  linear  analysis  techniques. 

Y-l  curves  for  each  of  the  nonlinear'  deviees  are  assumed  to  be  available,  so  for  any  one 
vector  v  the  vector  £(v)  may  be  evaluated.  This  may  require  some  interpolation  if  three  or 
more  terminal  nonlinear  devices  are  present.  Front  a  knowledge  of  v  and  £(v),  by  use  of  Kq.  (1) 
the  corresponding  vector  V  can  be  found.  Also  for  any  one  vector  v  a  unique  vector  Yg  exists. 
The  reverse  is  also  true.  Thus  given  sufficient  space  a  large  enough  table  can  be  constructed 
that  will  permit  reverse  interpolations,  i.e.,  given  a  vector  Yg  the  corresponding  vector  v  can 
be  found  to  within  some  specified  accuracy. 

Now  consider  Fig.  2(a).  For  any  value  of  v,  v.  a  corresponding  unique  value  f(v.)  exists. 
Since  f(v.)  is  a  current,  a  resistance  R.  v./f(v.)  may  be  defined.  Again  for  a  given  v.,  R.  is 
unique.  So  f rom  a  given  Y-I  curve,  R.  can  be  easily  found  for  each  v.  The  same  is  quite  obvi¬ 
ously  true  of  a  family  of  curves  [Fig.  2(b)];  however,  here  two  variables  v.  and  v  .  are  required 
to  define  R.. 

So  for'  any  given  vector,  V  ,  each  nonlinear  element  can  he  replaced  hy  its  corresponding 
value,  R..  In  effect,  for  each  value  of  Vg  the  entire  network  can  be  replaced  by  a  purely  positive 
linear  network.  If  we  consider  the  2 -port  N  (Fig.  1)  and  assume  that  all  its  Y  .  are  on  ports 
other  than  1  and  2,  then  for  each  value  of  Y  one  can  write  down  any  of  the  2 -port  matrices  by 
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straightforward  linear  network  analysis  methods.  Consider  the  h-matrix;  then 


V1  hii(V 

+  h12(~s>  V2 

h21(Xs> 

M  +  h22(L>  V2 

(2 ) 

as  V  is  changed  to  some  new  value,  the  h-parameters  will  also  change.  This  dependence  is 
implied  bv  the  notation  used  in  Eq.(2). 

So  far  it  has  been  assumed  that  the  elements  of  Vg  are  DC  voltages.  But  the  main  aim  of 
this  paper  is  to  investigate  the  frequency  performance  of  the  networks  considered.  The  fore¬ 
going  discussion  can  be  applied  to  this  aim  in  the  following  manner.  Assume  that  eac  h  of  the 
voltages  in  Vg  is  some  sinusoidally  varying  signal.  The  frequency  of  each  of  these  voltages  used 
need  not  be  the  same.  It  will,  however,  be  assumed  that  the  magnitude  of  one  V  .  is  far  larger 
than  that  of  all  the  others.  The  reason  for  this  assumption  will  become  apparent  later,  l.et  this 
dominant  voltage  be  denoted  by  V  and  its  frequency  by  f  .  Now  let  V  and  all  the  other*  V 
be  sampled  over  one  complete  cycle  of  V  Each  set  of  samples  will  determine  one  complete 
vector  V  and  hence  one  value  for  each  of  the  h-parameters.  If  n  samples  are  taken,  this  will 
determine  the  n  different  values  each  of  the  h-parameters  assumes  over*  one  evcle  of  \  .  Since 

Y  was  assumed  to  be  dominant,  this  pattern  will  be  repeated  to  a  high  degree  of  accuracy  over* 
each  cycle  of  Ygo.  This  scheme  is  illustrated  in  Fig.  5,  taking  h^Yg)  as  an  example.  Of  course 

all  the  other  \  present  are  also  sampled  at  the  same  rate  as  V  .  The  envelope  of  the*  result - 
si  M  M  so  1 

ant  magnitudes  of  h..(V  ),  (i,  j  1,2)  represents  the  variation  of  these  parameters  in  time  over* 

1  j  s 


Pft-  6  t?r?s] 


(b) 


Fig.  5.  (a)  \  input  waveform,  (b)  possible  resultant  h^(Yg)  wavcfor  i 


a  period  T  .  The  reason  for  choosing  one  dominant  \  .  is  now  apparent.  If  more  than  one  \ 

r  o  H  si  si 

were  dominant  the  periodicity  of  the  h-parameters  could  well  become  indeterminate,  at  least  in 

a  general  case.  In  many  practical  networks  the  assumption  that  one  drive  is  dominant  as  far  as 

the  nonlinearities  are  concerned  is  founded.  The  example  to  be  considered  will  illustrate  this. 

At  this  stage,  cases  have  to  be  considered  on  an  individual  basis,  especially  if  the  dominance  of 

onlv  one  \  .  can  no  longer  be  assumed.  In  this  case,  each  of  the  four  h.  parameters  can  be  e\- 

si  *  ij  r 


panded  in  a  Fourier  series  with  a  fundamental  period 


Thus 


4 


on 


ttn  ,  exp|jiwot| 


(3a) 


hll(t)  I 

i  0 


h12(t)"  I  H12,i 


i  0 


h21(t)  -  H2  exp[jiu.ot] 


i  0 


h22("  L,  H22,i 
i  0 


(  3h) 
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So  the  problem  posed  at  the  beginning  has  now  been  reduced  to  a  2 -port  problem  with  a  set  of 

h -parameters  given  in  ICq.  (3).  The  amplitudes  and  fundamental  frequency  of  these  parameters 

is  determined  b\  \  A  number  of  other*  Y  mav  also  be  present.  It  has  to  be  noted  that  an\ 

so  si  r 

DC  voltage  supplies  can  also  be  represented  by  a  \  . .  These  DC'  voltages  may  be  any  required 

magnitude,  comparable  or  even  greater  than  \  ,  since  they  will  not  introduce  a  periodic 

variat  ion. 


Kig.  6.  The  2 -port  N  shown  in  Fig.  1 
with  source  and  load  terminat ions. 
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In  the  most  general  case,  the  2-port  considered  will  be  supplied  by  some  voltage  \  ai  fre  - 
quonc\  f  (Fig.  b)  on  port  1  and  work  into  some  load  .  The  voltages  and  <  urrents  on  the  two 
sides  of  the  network  will  have  frequency  components  at  (f  ±  kf  )  for  k^  so  denoting 

v2i.  csPl.i  by  \'2(.  and  |ilm|  exp|j  by  llm, 

0© 

v2(U  ^  V  exP Life's  + 

(  -  on 

on 

i ,  (t )  v  1.  exp  [  \(jl'  +  mw  )  t|  .  Mb) 

1  —i  lm  1  IJ  s  o 

m  -  °° 

Similar  expressions  holt!  for  the  other  parameters.  Thus  from  Kqs.  (2),  (3)  and  (4): 
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+  ^  11.,  .  expfiio;  t]  x  N  \  exp[ju;  4  iu ;  ) 

^  12,1  ‘  oJ  Lj  2c  r  1  •'  s  o 
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x  V  I  S' 
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but  since  v^(t)  will  also  be  of  the  same  form  as  the  expressions  in  Kq.(4),  r.e., 


(5) 


V1U'  '  S  X  In  cxPl.'<“s  +  n“0)  ll 
n  -« 

it  follows  from  Eqs .  (8)  and  (6)  that 
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V,  Y  [II,,  .(I,  +1,  .)  +  H,,  (\  ,  .  +  V,  ,  .)| 

In  -j  1  ll,i  l,n-t  l,n  +  i  12, l  2,n-i  2,nti  1 
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and  by  the  same  argument 
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*2 1  ,  i^l  ,n-i 


+  1 


1  ,n4i 


)  4-  \\ 
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Since  — 00  <  n^  «  and  — 00  k  <r  <*,  Eqs. (7)  and  (8)  define  a  2n  x  2k  matrix  with  both  dimen¬ 
sions  stretching  from  — 00  to  400,  [n  every  practical  case  the  coefficients  1 1  ^  .  will  decrease  as 
i  increases,  also  [  Eq.  (4)]  the  coefficient  \  ^  and  I^m  will  decrease  with  an  increase  in  k,  so 
only  a  finite  portion  of  the  matrix  lias  to  be  considered  for  any  required  accuracy.  Eor  example, 
assume  that  only  the  first  three  H-parameter  coefficients  need  be  considered,  then  Kq.(7)  re¬ 
duces  to 


\  ,  II.,  1,  HI..  .(I.  ,4  1,  , )  4  H,,  l\,  _  4  I, 

l,n  ll,o  l,n  11,1  l,n-l  l,n4l  11,2  l,n-2  l,n42 

+  II,  ,  ,(  l,  _  4  I ,  -  )  +  H  ,  —  \  4H,-,,(\~  ,  4  \  _  ,) 

11,3  l,n-3  l,n43  12, o  2,n  12,1  2,n-l  2,n*1 


12,2  2,n-2  2,n42  12,3  2,n-3  2,n42 


(9) 


Equation  (8)  will  reduce  to  an  analogous  expression.  If  it  is  further  assumed  that  only  \ 
\  ^  r  •  \  1  2  are  interest,  then  Eq.  (9)  may  be  rewritten  in  matrix  form  as 
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(10) 


Since  the  number  of  voltages  and  currents  must  be  the  same,  only  the  portions  of  the  ma¬ 
trices  and  column  vectors  inside  the  dashed  lines  need  be  considered.  Equation  (8)  ma\ ,  of 
course,  also  be  written  in  a  matrix  form  analogous  to  Eq.(10).  Denoting  the  left-hand  side  of 
Eq,  (10)  by  V ^  the  two  truncated  column  vectors  on  the  right  side  by  n  and  \  ^  n»  and  using 
a  similar  notation  for  the  other  relevant  entities,  Eqs.  (7)  and  (8)  may  be  rewritten  as 


- 1 
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3 

_ J 

U12" 

‘li.n 

~2,n 

^21  |  ~22 

~2,n 

Each  of  the  H-parameters  is  an  n  x  n  matrix  of  constants,  thus  Eq.(ll)  describes  the 
h-matrix  of  a  linear  2  port  with  n  input  and  n  output  ports,  however,  since  each  \  ^  n>  etc.,  is 
associated  with  one  specific  frequency,  this  2  port  differs  fundamentally  from  a  conventional  lin¬ 
ear  network  in  that  each  port  is  associated  with  a  distinct  frequency  different  from  all  the  others. 
Figure  7  is  the  pictorial  representation  of  Eq.  (11).  Despite  the  fact  that  each  port  is  associated 
with  a  different  frequency,  the  network  itself  is  linear  and  hence  any  transfer  parameter  can  be 
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Fig.  7.  Resultant  linearization 
of  the  nonlinear  2-port  N. 


Fig.  8.  Balanced  mixer  terminated  in  matching  networks  and  N^. 
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computed  using  simply  linear  methods.  Thus,  for  example,  if  the  power  transfer  between  inpul 
port  1  and  output  port  i  is  required,  this  will  he  given  by  the  |S^  .  |  wher  e  .  is  the  1,  i  en¬ 
try  in  the  scattering  matrix  of  network  N^.  As  far  as  the  original  2-port  N  is  concerned,  |S^ 
defines  the  power  transfer  from  port  1  at  frequency  f^  .  to  port  2  at  frequency  f^  ^  Any  other 
parameters  may  be  evaluated  in  a  similar  manner. 

To  recap,  given  a  2-port  N  containing  nonlinear  resrstors,  frnd  the  time  performance  of  its 
h -parameters  (or  any  other  set  for  that  matter),  expand  these  in  a  Fourier  series,  collect  the 
dominant  coefficients  to  form  Fq.  (11)  and  the  nonlinear  2-port  problem  is  reduced  to  an  equiv¬ 
alent  linear  2n-port  problem.  For  example; 


The  Balanced  Ring  Mixer 

Problem  Investigate  the  frequency  performance  of  the  balanced  mixer*  (Fig.  8)  and  deter¬ 
mine  the  cha racteristics  of  matching  networks  placed  on  the  input  and  output  ports  (networks 
and  N^),  which  will  give  the  smallest  insertion  loss  between  the  signals  at  source  frequency 
and  the  required  output  frequency. 

be t  the  R F  signal  at  angular  frequency  be  applied  to  port  1,  the  LO  signal  to  port  3,  and 
let  port  2  be  the  output  port  at  which  the  IF  signal  at  angular  frequency'  (a;g  —  u.*^)  is  extracted. 

The  3-coil  t ransformers  are  assumed  to  be  frequency'  invariant  over  the  band  of  interest.  Assum¬ 
ing  the  FO  drive  is  dominant  (this  is  a  very  fair  assumption  in  the  case)  it  will  be  found  on  in¬ 
spection  that 


v .  v ,  V 
1  2  so 


v_  v  ■  —  V 
3  4  so 
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This  simplifies  matters  very  much  since  none  of  the  interpolation  procedures  described  at 
the  beginning  need  be  used.  It  is  also  assumed  that  the  diodes  are  matched,  and  hence,  have 
essentially  the  same  \  -1  characteristics.  These  will  be  taken  to  be  of  the  form  shown  in  Fig.^(ai. 
With  this  assumption  and  the  V  drive  as  shown  in  Fig.  9(h),  the  resultant  diode  resistance  per¬ 
formance  against  time  will  very  nearly  be  that  shown  in  Fig.  9(c).  \n  additional  degree  of  sophis¬ 
tication  is  added  in  Fig.  9(c)  by  assuming  the  rise  and  fall  times  are  unequal.  From  Kq.  (12)  it 
follows  that  diodes  1)^  and  1),  will  have  a  resistance  waveform  of  the  type  shown  in  Fig.  9(c). 

I  .et  this  resistance  be  denoted  by  R(t).  The  resistances  of  diodes  1)^  and  1)^  will  have  the  same 
shape  as  R(t),  however,  they  will  be  shifted  by  180°  with  respect  to  R(t).  For  this  reason  these 
impedances  will  be  referred  to  as  R  (t).  Bv  a  straightforward  Fourier  analysis  of  Fig.  9(c),  it 
follows  that 


where 


H+(t ) 

l(I?b 

+  Rr)  +  j(R,)  -  R,.)  < (t) 

R  _(t) 

bn  b 

! 

^•|rq 

1 

+ 

c(t>  “ 

or> 

v  1 

r  ,  fia  . 

sin  nrt  —  ) 

sin  nr(  t|t-  ) 

nr 
sin  — y 

x  ^ 

Lt  2 
n  1 

L  n7r(  — ) 

,  ah . 

nr(  tjt-  )  J 

nr 

2 

.  nr 
sin  — 

V  k(n)  - t~-  cos  nu;  t 

nr/2  o 

n  1 


(  Ha ) 


( 1  3b) 


(  14) 


(14a) 
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I  =  f (v) 


RESISTANCE 


I II- 6- 12128 1 


Fig.  9.  (a)  Diode  characteristic,  (b)  Vso  cosine  wave  versus  time,  (c)  resultant 

resistance/time  performance  of  a  diode  driven  by  the  Vgo  shown  in  (b). 


This  is  a  trapezoidal  wave  of  ±\  magnitude.  It  should  be  noted  that  if  <5a  6b  0,  k(n)  1,  then 
for  all  n  and  t(t)  reduces  to 


u(t)  -  ^  - —  cos  no?  t  (15) 

— i  nr  o 

n  1  2 

which  is  the*  Fourier  expansion  of  a  square  wave  varying  between  +1  and  -1. 

An  equivalent  derivation  on  a  conductance  basis,  assuming  the  two  diode  slopes  to  be  Gj  and 
Cif  gives 


G  (t)  - 

l(Gf  +  C.b) 

1 

M  ^ 

r\ 

~r 

1 

c(t) 

(16a) 

Gjn 

2  (°f  *  (:h’ 

1  ?‘Gf-Gb» 

€(t)  . 

(  IC>b) 

G  +  (t)  and  G  (t)  are  defined  in  the  same  fashion  as  H^(t)  and  R  (t). 
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Assuming  the  instantaneous  resistances  (conductances)  of  the  four  diodes  in  Fig.  8  to  be, 
respectively,  R^,  R^,  R^,  R^  (G^,  G^,  G^,  G^)  then  the  instanteous  h-matrix  between  ports  1  and 
2  is  given  by 


(h) 


hf1  h12 


h21  h22 


( - j -  +  - 5 - ) 

'Gd  f  (,3  G2+G4' 


R1R2~R3R4 

(K,  t  R3)  (R2  +  R4> 


RIR2 


R3R4 


(Rd  +  Rj) (H  +  R4I 


(k4  +  «3  +  R2  *  0 


(171 


but  R  R.  -  R  (t)  and  IF  =  R„  ~  R  (t),  or  vice  versa. 
12+  3  4- 


A  similar  set  of  relations  holds  for  the  Y  parameters.  Using  Fqs.  (13),  (16)  and  (17). 
2 


h 


1 1 


°b  4  Gf 


( 1  8a  I 


12 


22 


-h 


21 

2 


Rb  +  Rf 


(  I  8b  I 


(  18cl 


It  is  interesting  to  note,  that  and  h  reduce  to  time  invariant  immittanres.  I  sing  <  (i  ) 


in  the  form  given  in  Kq  (14a) 

M 


12 


21 


R.  -  R,  \  f  a 

TT,^  TTt~)  [k(1)  cos  wo‘  -  T  k(3)  003  3<*V  +- 


(2n  V  1 )  k(2n  +  i]  cos(2n  f  i]  woG 


follows  that 


H 


11,0 

Gb +  Gf 

H,  ,  ■ 
11,1 

i^O 

2 

1L, 

22,1 

22,0  = 

Rb  +  Rf  ’ 

1^=0 

12,o 

"21.0  0 

=  0 


=  0 


(  1  Pa  > 

(  IPbi 
(  1  9c  i 


II 


12,  i 


i^o 


=  II 


21,  i 


i960 


2  /Rb  Hf\  v 

»  '  Rb  +  Hf  ' 

1=1 


SitlT 


k(i) 


( 1  9d ) 


liquation  (19)  completely  defines  all  parameters  in  the  2n-port  H-matrix.  Table  1  gives  a 
number  of  ratios  for*  k(n).  In  each  instance  the'  rise  and  tall  times  have  been  assumed  equal 
(i.e.,  <5^).  It  should  be  noted  that  the  fundamental  component  (i.e.,  n  1  )  in  each  of  the  ease 

considered  is  very  close  to  unity,  and  as  n  increases  k(n)  departs  from  unity.  So  the  higher 
order  harmonics  arc  attenrrated  to  a  larger  degree  than  the  fundamental.  This  fact  will  actually 
aid  the  mixer  design.  More  will  be  said  about  this  when  the  numerical  results  of  thrs  example 
are  diserrssed. 
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TABLE  I 

RATIOS  FOR  k(n) 

6  -  6.  -  0.01 

6  =  6.  =  0.05 

6  6,  -  0.1 

n 

a  b 

a  b 

a  b 

1 

0.999812 

0.995865 

0.983641 

3 

0.998433 

0.963397 

0.858397 

5 

0.995865 

0.900316 

0.636620 

7 

0.991945 

0.810335 

0.367884 

9 

0.986725 

0.698648 

0. 109294 

Assuming  6^ 


0,  the  H-matrix  of  the  same  order  as  that  in  Kq.  (10)  becomes 
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22.0.0 
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0 
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1 
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0 
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(20) 


where  II.  .  ,  1  (Rt  +  RJ  evaluated  at  cc  -  w  +  ka;  and  .  l/(H,  +  KJ  evaluated  at 

1 1 ,  o ,  k  of  so  22,0, k  b  f 

a’  <jc s  +  ku?^.  Equation  (20)  represents  a  10-port  passive  reciprocal  network  with  5  input  and 
S  output  ports.  A  number  of  important  conclusions  about  the  character  of  this  10-port  can  be 
drawn  from  Kq.  (20).  This  can  best  be  done  by  writing  out  some  of  the  equations,  thus: 


1,-2 

nH  ,  i-  , 

1  l,c>,  —2  1,-2 

+  '‘  2,-1  _  T  x  2,1 

(21a) 

1.-1 

IK  4  aK  A 

1 1 ,o, -1  1,-1 

+  X  2,-2  +  X  2,o  “ 

T  V2,2 

(211)) 

Vl,o 

IK,  I,  + 

11,0,0  1,0 

X  2,-1  +  X  2,1 

(21c) 

V1.1 

"l  1,0,1  1 1  .  1 

T  V2,-2  +  X  2,o  + 

X  2,2 

(21d) 

V1,2 

II, ,  ,  1,  ,  - 

1 1,0,2  1,2 

TV2,-1+X2,1 

(21c) 
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The  frequency  associated  with  any  one  port  is  determined  by  the  voltage  (or  current  suffix), 
thus  \  ^  .  is  at  a  frequency  u,'g  +  itr^,  etc.  So  from  Eq.  (21)  it  follows  that  components  at  even 
multiples  of  the  angular  frequency  cj ^  on  the  output  ports  will  produce  components  at  odd  mul¬ 
tiples  of  u;q  at  the  input  ports  and  vice  versa.  From  1  ig.  8  it  can  be  seen  on  inspection  that  the 
input  voltage  components  will  contain  the  voltage  ()  (i.e.,  at  u;g),  so  there  will  be  a  current 
component  1^  .  This  will  result  in  only  odd  suffix  voltages  and  currents  on  the  output  pons, 

which  in  turn  will  produce  only  even  suffix  voltages  and  currents  on  the  input  ports.  As  a  result, 
only  even  harmonies  of  lo  will  be  present  on  the  input  ports  and  only  odd  harmonic's  on  the  out¬ 
put  por  ts.  This  has  an  impor  tant  consequence  because  it  implies  that  all  .  for'  r  odd. 

and  V ^  ^  \  ^or  *  even,  will  have  zero  amplitude.  So  two  of  the  input  ports  and  three  of  the 

output  ports  of  the  10  ports  are,  for  practical  purposes,  nonexistent.  So  1  hey  may  be  ignored. 
Since  the  mixer'  will  be  used  (in  this  example)  as  a  down  converter,  the  component  \  ?  ^  is  the 

required  output.  So  from  equations  of  the  type  given  in  Eq.  (21)  the  input  and  output  voltages,  in 
decreasing  order  of  magnitude,  can  be  written  as  follows: 


Input  Magnitudes 


Output  Magnitudes 


1  ,o 


2,-1 


1,-2 


2,1 


1,2 


2,-3 


1,-4 


2,3 


1,4 


2,-5 


and  similarly  for  currents.  I’sing  these  magnitudes  as  the  new  column  vectors,  a  new  10-pori 
matrix  may  be  constructed  by  inspection  of  Eq.  (10),  thus 
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By  inspection  of  Fig.  8  it  follows  that  Fq.(lZ),  which  governs  the  voltage  distribution  across 
the  diodes,  will  be  unaffected  by  connecting  general  2-ports  ( N ^  and  N  )  in  series  with  pons  1 
and  2.  If  this  is  done,  the  10-port  can  be  represented  as  shown  in  Fig.  10. 


j;s  rum 

INPUT  PORTS  OUTPUT  PORTS 


Fig.  10.  Augumented  10-port  network. 

The  networks  and  appear-  on  each  port  evaluated  at  the  port  frequency.  Let  the  aug¬ 
mented  10-port  input  magnitudes  be  denoted  by  and  ^  ^  and  those  on  the  output  ports  by  \  j 
and  (Fig.  10),  then  since  the  transmission  matrix  of  is 


the  overall  transmission  matrix  of  networks  N.  .(  I  .)  becomes 

l,i  ~  1 


V  * 

A,  '  B, 

V  ,  . 

~  A 

~  1  i  ~  1 

__L_ 

— 1,1 

I  A 

c\  !  i). 

1 , 

-A 

M  i  -i 

L  i  J 

1 , 1 

where  each  of  matrices  A.,  B.,  etc.,  is  a  diagonal  mat  i  ix  with  entries  A.  ,  \,  etc.,  for 

11  H  l.o  1,-2 

( A ^ )  in  the  leading  diagonal.  A  similar  matrix  (T^)  may  be  defined  for  the  output  ports,  thus 
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(2*4) 


\\  . 

A-  j  B_ 

-2,1 

—2  ~2 

_ i _ 

-B 

-2,  i 

i7\i"z 

L  i  J 

-B 

Converting  the  11-matrix  of  the  unaugmented  10-port  into  an  equivalent  transmission  matrix 
(T  ),  then  the  overall  transmission  matrix  of  the  augmented  10 -port  (T)  is  given  by 


|T|  [T.J  [T  |  [TJ 


(251 


If  it  is  now  assumed  that  the  load  and  source  resistor's  R  and  R.  are  equal,  and  1  hat  R 

s  \j  s 

Rj  1,  then  the  scattering  matrix  parameters  of  the  augmented  network  is  given  b\ 

-1 


S14  U  -  2  x  (C  +  D)  x  A 


Ubn  > 


-1 


^2^21  -2XA 


s  22  r  -  2  x  a'1  x  <b  +  o) 


(261)i 
(26c  i 


has  the  property  that 
5 16 


unit  matrix  and  A 

A  t  B  +  C  +  1).  The  scattering 

mat  rix 
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Max .  available  power*  at  input  port  1 

Power-  delivered  to  frequency  (u,*g  —  u.*  ) 
.Max.  available  power  at  source  frequency  **.• 


(27  1 


(28» 


In  similar*  fashion  all  the  other  j  s  ^  ^  arr  r*atios  of  the  powers  delivered  to  the*  frequency  ai 

por-t  j  to  the  maximum  available  power*  from  the  source.  So  an  evaluation  of  the  permils  a 
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computation  of  how  much  power  goes  into  each  of  the  frequencies  present  on  both  sides  of  the 


mixer. 


Also,  the  input  impedance  Z^  at  the  j 


.th 


port: 


Z.  = 
J 


1  —  s. 


r~r 


jj 


(29) 


So  from  a  knowledge  of  the  s^,  the  input  impedances  on  both  sides  of  the  mixer  at  all  the  fre¬ 
quencies  are  available. 

Using  the  above  mathematical  background,  a  computer  program  has  been  written  to  evaluate 
these  power  and  impedance  distributions. 


MODCAP 

Based  on  the  mathematical  background  described,  a  Fortran  program  was  written  to  analyze 
the  performance  of  a  modulator  connected  between  two  general  filters,  hence  the  name  MODCAP 
for  MODulator  Circuit  Analysis  Program.  The  MODCAP  program  appears  at  the  back  of  this 
report . 

The  purpose  of  the  program  is  to  evaluate  the  power  transfer  ratio  from  the  input  port  at 
input  angular  frequency  to  the  power  at  all  the  other  frequencies  assumed  present  at  both 
ports.  For  example,  if  the  strongest  five  harmonics  at  the  input  and  output  ports  arc  considered, 
the  situation  depicted  in  Fig.  10  ensues.  The  program  will  evaluate  the  power  transfer  (as  a  ratio 
in  dB)  between  input  port  1  and  any  of  the  remaining  nine  ports.  It  will  also  evaluate  the  input 
impedances  at  the  input  and  output  ports  at  all  the  frequencies  considered.  Further,  the  program 
will  evaluate  the  amount  of  change  induced  in  the  power  transfer  ratios  and  impedances  if  some 
element  or  elements  constituting  the  network  change  by  a  specified  amount.  All  of  these  cal¬ 
culations  can  be  carried  out  in  three  different  modes: 

1 .  A  single  frequency  is  specified  in  the  data  for  the  input  and  local  oscillator 
drive. 

2.  A  band  of  input  frequencies  around  some  nominal  value  is  specified  together 
with  a  single  local  oscillator  frequency.  In  this  case,  the  output  frequencies 
are  also  grouped  in  bands. 

3.  The  input  frequency  extends  over  a  band  and  a  single  local  oscillator  fre¬ 
quency  is  chosen  to  be  a  single  frequency.  The  program  will  automatically 
change  the  nominal  local  oscillator  frequency  for  each  frequency  in  the 
input  band  such  that  a  specified  single  output  frequency  results. 

In  each  case  the  program  evaluates  the  harmonics  automatically  at  single  frequencies  or  in 
frequency  bands  depending  on  the  mode  used.  In  its  present  form,  due  to  storage  limitations,  the 
program  has  these  maximum  values: 

1 .  The  filters  terminating  the  mixer  on  each  side  must  have  no  more  than 
10  nodes  and  contain  only  two  or  three  terminal  elements.  For  ten  nodes 
there  can  only  be  three  terminal  elements. 

2.  The  maximum  number  of  harmonics  on  each  side  of  the  mixer  cannot 
exceed  10. 

3.  The  same  number  of  harmonics  always  have  to  be  considered  on  each  of 
the  two  ports . 

4.  If  frequency  bands  around  each  harmonic  are  used,  they  can  contain  no  more 
than  21  frequencies. 

5.  If  effects  of  one  or  more  clement  changes  arc  investigated,  only  42  sets  of 
variations  can  be  accommodated. 

All  of  these  constraints  can  be  overridden  by  appropriate  dimensioning. 
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Data  is  entered  via  a  subroutine.  Data  for  a  typical  case  are  shown  in  the  MODCAP  input 
data  (program)  at  the  end  of  this  report  where: 

M  specifies  the  number  of  frequencies  in  each  harmonic  band. 

N  is  the  number  of  harmonics  including  the  fundamental  considered  on  each  side. 

H1NOD  (a)  is  where  a  =  1  and  2  for  the  highest  node  numbers  in  network  1  and  2, 
respectively. 

RISE  is  the  fractional  rise  time  of  switching  wave  shape  (6^/T)  (Fig.  9c). 

FALL  is  the  fractional  fall  time  of  switching  wave  shape  (6^/T)  (Fig.  9e). 

RT  is  the  magnitude  of  source  and  load  resistor.  Tf  unequal  terminations  arc 
required  transformers  have  to  be  included  in  networks  1  or  2. 

X(a,  i,  j)  is  element  X(L,  R  or  C)  of  network  a(l  or  2)  appearing  between  nodes 
i  and  j.  The  units  are  Henries,  Ohms  and  Farads  and  the  values  are  entered 
in  floating  point  form,  e.g.,  L(l,  3,  5)  =  5.86E-6  in  network  1.  The  inductor 
appearing  between  nodes  3  and  5  has  a  value  of  5.86  pH. 

Tf  3- terminal  devices  (e.g.,  transformers)  appear  in  either  of  the  terminating  networks, 
their  presence  and  location  is  described  as  follows: 

Each  device  is  assigned  a  consecutive  number.  Their  input  node  is  labelled  El, 
the  output  node  E2  and  the  common  node  E3.  The  network  number  and  device 
are  in  brackets: 

El (2  1 )  -  1 
E2(2 ,  1  )  -  2 
E3(2 ,  1  )  -  3 

This  notation  tells  the  computer  that  in  network  2  the  input  node  of  3- terminal  device  1  is 
connected  to  the  node  labelled  1.  The  output  node  of  device  1  in  network  2  is  connected  to  node  2 
and  the  common  node  of  device  1  in  network  2  is  connected  to  node  4.  This  describes  the  loca¬ 
tion  of  device  1  in  network  The  electrical  parameters  of  3- terminal  devices  are  specified  by 
their  H- matrix.  All  four  H- parameters  must  be  frequency  independent,  fill  and  H22  are  in 
units  of  Ohms  and  Mhos,  respectively;  the  remaining  two  are  dimensiomess .  All  magnitudes 
are  again  specified  in  floating  point  form.  For  example,  HI  1(2,  1)  =  0.5EO  means  that  Hll  of 
device  1  in  network  2  is  0.5  Ohms. 

If  frequency  bands  are  used  these  arc  entered  as  Fll(i)  -  value  where  i  1,  2,  ...  M.  If  the 
frequencies  are  equally  spaced  they  can  be  entered  by  specifying  the  first  frequency  KR(1) 
'value',  the  increment  DR(2)  -  'value'  and  the  last  frequency  FR(7)  =  'value'.  Frequencies  can 
be  entered  sequentially  or  with  increments  or  a  combination  of  both.  There  mav  be  as  many 
increments  as  necessary,  but  the  total  number  of  frequencies  cannot  exceed  21. 

FLON  is  the  nominal  local  oscillator  frequency. 

FIFN  is  the  nominal  output  frequency.  This  is  0.0  for  modes  1  and  2  and  assumes 
some  non-zero  value  only  for  mode  3. 

This  completes  the  nominal  data  input.  If  no  variations  are  required,  a  RETURN  card 
follows  FIFN  and  then  'i  CONTINUE'  cards  where  i  -  1,  2.  .  .  42.  After  42  CONTINUE,  a  STOP 
and  an  END  card  completes  the  data  input. 

If  variations  are  required,  the  2  CONTINUE  card  follows  the  RETURN  card.  A  'VAR  1' 
card  is  inserted  next.  Following  is  the  first  set  of  variations.  Thus  NEWR(1,  1,5)  =  0.607 EO 
means  that  the  old  value  of  the  resistor  in  network  1  between  nodes  1  and  5  assumes  the  new 
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value  of  0.607  Ohms.  At  the  end  of  the  first  set  of  variations  a  RETURN  card  has  to  be  inserted. 
If  there  is  more  than  one  variation,  the  procedure  is  repeated  with  VAR  =  2,  etc. 

NOTE: 

The  nodes  in  the  two  terminating  networks  have  to  be  labelled  as  follows 
the  input  node  has  to  be  node  4,  the  output  node  is  labelled  2,  internal  nodes 
are  labelled  arbitrarily,  and  the  highest  node  number  is  assigned  to  the 
common  ground  node. 
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